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ABSTRACT 

Using the Atacama Large Millimcter/submillimcter Array (ALMA), we have conducted a blind 
rcdshift survey in the 3 mm atmospheric transmission window for 26 strongly lensed dusty star- 
forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected 
to have Si.4 mm >20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio 
(S843MHz < 6mJy) or far- infrared counterparts (Sioo^m < 1 Jy, S60/im < 200 mJy). We robustly detect 
44 line features in our survey, which we identify as redshifted emission lines of 12 CO, 13 CO, C i, H2O, 
and H2 + . We find one or more spectral features in 23 sources yielding a ^90% detection rate for 
this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single 
line. For the sources with only one detected line, we break the redshift degeneracy with additional 
spectroscopic observations if available, or infer the most likely line identification based on photometric 
data. This yields secure redshifts for ^70% of the sample. The three sources with no lines detected 
are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our 
sample is z=3.5. This finding is in contrast to the rcdshift distribution of radio-identified DSFGs, 
which have a significantly lower mean redshift of z=2.3 and for which only 10-15% of the population 
is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution 
and compare our measured redshift distribution to that of models in the literature. 
Subject headings: cosmology: observations — cosmology: early universe — galaxies: high-rcdshift 
galaxies: evolution — ISM: molecules 
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1. INTRODUCTION 

In the last decade, impressive progress has been made 
in our understanding of galaxy formation and evolution 
based on multi-wavelength deep held studies. Millimeter 
(mm) and submillimctcr (submm) continuum observa- 
tions demonstrated that luminous, dusty galaxies were 
a thousand times more abundant in the early Universe 
than they are a t present day (e . g.. iSmail et al.l 119971 : 
iBlain et all 119991 : IChapman eHd] 120051 ). The hrst sur- 
veys of the redshift distribution of dusty star-forming 
galaxies (DSFGs) suggeste d that the DSFG populatio n 
peaks at redshift ~ 2 (e.g.. IChapman et al.ll2003l I2005D . 
coeval with the peak of black hole accretion and the peak 
of the star for mation rate density as mea sured in the op- 
tical/UV (e.g.. lHopkms fc Beaconll2006l ). These studies 
suggested that the bulk of star formation activity in the 
universe at z = 2 — 3 could be taking place in DSFGs, 
hidden from the view of optical/U V observations due 
to the high dust obscuration (e.g.. iHughes et all 119981 : 
IBlain et al.l [19991) . 

Optical surveys now allow estimates of the history of 
star formation ( the ' Madau-Lillv' plot: iMadau et al\ 
119961: iLillv et alj|1996t iHopkins fc Beacoml 120061) out to 
z ~ 8 (e.g.. iBouwens et al.l 120101 l201uV but have un- 
certain dust extinction corrections. Submm observations 
can provide a more complete picture of the amount of 
highly obscured star formation over a large range of 
look-back times. However, such studies have been ham- 
pered by the difficulty of obtaining robust redshifts for 
DSFGs. This difficulty increases strongly as a function 
of redshift, and mainly arises from the coarse spatial res- 
olution (-~20") of single-dish submm observations and 
the dust-obscured nature of the sources, which often 
prohibits identification of counterparts at other wave- 
lengths. The solution has been to obtain higher spa- 
tial resolution data, usually at radio and/or mid- infrared 
wavelengths, in which the most likely count erpart to the 
submm emission could be identifi ed (e. g. 

Ashbv et al.l[ 2006: Po pe et al.ll20olTWardlow 



2002; 



2011 



et al. | 
ct al. 



Yun et al.ll2012l ). The slope of the spectral energ; 



distributions (SEDs) of galaxies in the radio or mid- 
infrared (MIR), however, is such that the K-correction 
is positive, and galaxies become more difficult to detect 
at high redshifts. By contrast, the steeply rising spec- 
trum of dusty sources leads to a negative K-correction 
for DSFGs at submm wavelen gths, resulting in fluxe s 
roughly constant with redshift ( Blain fc Longairlll993l ). 
Therefore, while DSFGs may be discoverable at submm 
wavelengths at almost any redshift, their emission may 
be hidden at other other wavelengths. Indeed, in submm 
surv eys typically 50% o f DSFGs lack robust counterparts 
(e.g. jBiggs et al.ll2011l ) albeit the fraction depends on the 
depth of the radio/MIR observations. This mismatch in 
the wavelength sensitivity could bias the DSFG redshift 
distribution, particularly at z > 3. 
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A more reliable and complete method to obtain se- 
cure multi- wavelength identifications of DSFGs is to fol- 
low the single-dish detections up with mm intcrfcrom- 
etry. Prior to ALMA this method has proven to be 
time-intensive, requiring entire nights of time detect a 
single source; the first sample detected blindly in the 
continuum with mm interferometry was published by 
lYounger et al.l (120071). A l arger sample was published 
recently by iSmolcic et al] ([2012D . which included op- 
tical spectroscopic redshifts for roughly half the sam- 
ple and photometric redshift estimates the remaining 
sources in the sample which suggested that the previous 
spectroscopically dete rmined redshift distributions (e.g., 
IChapman et al.|[2005l) were biased low. 

A more direct and unbiased way to derive redshifts of 
DSFGs is via observations of molecular emission lines at 
millimeter wavelengths which can be related unambigu- 
ously to the (sub)mm continuum source. This method 
has only become competitive over the past years with 
the increased bandwidth of mm/submm facilities. Its 
power to measure reliable redshifts has been demon- 
strated in the case of SMMJ14009+02 52 and HDF850.1 
(jWeifi et al.ll2009al : IWalter et all 120121 ). two of the first 
DSFGs detected by SCUBA, for which other methods 
failed to deliver redshifts for more than a decade. While 
CO redshift surveys of a representative sample of DSFGs 
will remain observing time expensive till the operation 
of full ALMA, CO line redshifts f or strongly lensed sys- 
tems can be obtained easily (e.g. iSwinbank et al.1 120101 : 
ICox et al.ll20~L"Ll iFraver et al.1120111 ). 

In the past studies of strongly lensed sources have 
been limited to a handful of targets due to their 
rareness and the lack of large scale mm/submm sur- 
veys. This has changed dramatically over the past 
years with the advent of large area survey s from Her- 
schel ( specifically H- AT LAS and HerMES lEales et al.1 
[20101: lOliyer et al.|[20Tfj|) a nd the South Pole Telescope 
(SPT-SZ. lCarlstrom et al.ll2011l) . These surveys have de- 
tected hundreds of strongly lensed high -redshift DSFGs 
(IVieira et al.ll2010HNegrello et al.ll2010D Firs t CO red- 
shift mea surements at mm ( Lupu et al.1 120121 ) and cen- 
timeter ([Harris et al.l 120121 ) wavebands of H-ATLAS 
sources suggested that the lensed DSFGs lie within 
the sam e redshift range as unlensed, radio-identified 
sources ([Chapman et al.ll2005|) . Although a large over- 
lap between the SPT and Herschel populations is ex- 
pected, SPT's longer selection wavelength of 1.4 mm 
predicts a broader redshift distribution than Herschel 
detected sources and indeed photometric redshifts of 
DSFGs discovered by the SPT confirm this expectation 
([Greve et al.l[2012l) . 

In this paper, we present the results from an ALMA 
CO redshift survey of a sample of 26 strongly lensed 
DSFGs selected from 1300 deg 2 of SPT-SZ survey data 
(|Carlstrom et al.|[20TTh . The depth of the SPT-SZ survey 
data, which is sufficient to detect 5i.4 mm '^20 mJy sources 
at 5(7, combined with the flat redshif t selection func- 
tion o f DSFGs at this wavelength (e.g.. IBlain fc Longairl 
Il993f ). has produced an optimal sample for mm molec- 
ular line redshift searches in strongly le n sed D SFGs. 
In an accompanying paper, IVieira et all (I2013D show 
that these sources ar e virtually all strongly lensed, while 
iHezaveh et al.l (120131 ) report the associated lens modeling 
procedure. 
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In Section we describe the target selection and ob- 
servations. The biases on the observed redshift distribu- 
tion, resulting from the source selection and the effect 
of gravitational lcnsing are discussed in Section |4] Our 
results are summarized in Section [5] Throughout this 
paper we have adopted a flat WMAP7 cosmology, with 
H Q = 71 kms^Mp c" 1 , tt A = 0.73 and Q m = 0.27 
(jKomatsu et al.ll201lD . 

2. OBSERVATIONS 

We observed a sample of 26 bright (Si.4 mm > 20mJy), 
1.4 mm selected SPT sources with ALMA. Sources were 
selected from the first 1300 deg 2 of the now complete 
2500 de g 2 SPT-SZ survey (for more details on the sur- 
vey, see !MflhimioniFil][20il [Stor^iFin EoTJ) . The 
flux density cut is done on the initial raw flux den- 
sity, and not the final de-boos ted flux dens i ty, th e de- 
tails of which can be found in IVieira et al.l (]2010f ) and 
ICrawford et "all (l2010h . To remove synchrotron domi- 
nated systems we required dust- like spectra betwee n 1.4 
and 2 mm (Si.4 mm /S2.omm >2; IVieira et al.l [2010h . In 
addition, we used far-infrared (FIR) and/or radio cri- 
teria to remove low-redshift contaminants (see Section 
14.11) . In order to refine the relatively coarse SPT source 
positions (the SPT's beam size is 1'.05 at 1.4 mm) we 
further required follow up observations at higher spatial 
resolution (typically 870/^m images from the Large Apex 
BOlometer CAmera (LABOCA) or 1 mm data from the 
Submillimeter Array). Based on 1.4 mm flux densities, 
our Cycle targets comprise a representative sample of 
the SPT sources meeting these selection criteria. This is 
shown in Appendix ICl where we present the SPT 1.4 mm, 
LABOCA 870^m, and Herschel-SPIKE 350/Ltm flux den- 
sity properties of this subsample compared to all SPT 
sources which have been observed with Herschel and 
LABOCA. 

In order to optimize the ALMA observing efficiency, 
we assembled 5 groups of targets that lie within 15° 
of each other on the sky - this restriction precluded 
a complete flux-limited sample. We excluded two 
sources with redshifts previously determined by Z-Spcc 
(a wide-band, low resol ution spectrometer op erating be- 
tween 190-310 GHz, see Bradford et al.|[2004] ) on the At- 
acama Pathfin der Experiment (A PEX) telescope and 
XSHOOTER (jVernet et al.l 12011 o r the the FOcal 
Redu cer and Spectrograph (FORS2; lAppenzeller et al.l 
\l99B) on the ESO Very Large Telescope (VLT). 

The ALMA observations were carried out in 2011 
November and 2012 January in the Cycle early sci- 
ence compact array configuration. We performed a spec- 
tral scan in the 3 mm atmospheric transmission window 
with five tunings in dual polarization mode. Each tun- 
ing covers 7.5 GHz in two 3.75 GHz wide sidebands, 
each of which is covered by two 1.875 GHz spectral win- 
dows in the ALMA correlator. This setup spans 84.2 to 
114.9 GHz (with 96.2 to 102.8 GHz covered twice; see 
Figure [T]), nearly the entire bandwidth of the Band 3 
(84—116 GHz) receiver. Over this frequency range 
ALMA's primary beam is 61" — 45". The observations 
employed between 14 and 17 antennas in different ses- 
sions, and resulted in typical synthesized beams of 7" x 5" 
to 5" x 3" (FWHM) from the low- to high-frequency ends 
of the band. Each target was observed for ~ 120 seconds 
in each tuning, or roughly 10 minutes per source in total, 



not including overheads. 

Typical system temperatures for the observations were 
T sys = 60 K. Flux calibration was performed on planets 
(Mars, Uranus, or Neptune) or Jupiter's moons (Cal- 
listo or Ganymede), with passband and phase calibration 
determined from nearby quasars. The data were pro- 
cessed using the Com mon Astronomy Software Applica- 
tion package (CASA, iMcMullin et al.ll2007t iPetrv et al.l 
I2012D . Calibrated data cubes were constructed with 
a channel width of 19.5 MHz (~ 50 — 65 kms -1 for 
the highest and lowest observing frequency). The 
typical noise per channel is 2mJybeam _1 across the 
band and 1.4mJybeam _1 between 96.0 and 102.8 GHz 
where two tunings overlap. Continuum images gener- 
ated from the full bandwidth have typical noise levels of 
70 /iJy beam -1 . 

The spectral coverage of this experiment includes 
CO(1-0) for 0.003<2<0.36 and one or more CO 
lines, between the (2-1) and (7-6) transitions, be- 
tween 1.0<z<8.6, with the exception of a small redshift 
"desert" between 1.74<z<2.00 (see Figure[2]). An addi- 
tional redshift desert at 0.36<z<1.0 is also present, but 
our high 1.4 mm flux density threshold effectively re- 
quires that our sources be gravitationally lensed ( 34. 1|) 
and it is highly unlikely that sources at this redshift will 
be lensed {MM- 
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Fig. 1. — Spectral setup and frequency coverage of our 5 tun- 
ings (shown in different colors) in ALMA band 3 for the source 
SPT0103-45. In each tuning, four spectral windows covering 
1.875 GHz each were placed in contiguous pairs in the lower and 
upper sidebands (LSB/USB). Note that the frequency range 96.2- 
102.8 GHz (delimited by dotted vertical lines) is covered twice. The 
total spectral range is indicated by the dashed vertical lines. The 
top panel shows the atmospheric transmission across band 3 at 
Chajnantor for 3mm prccipitable water vapor (PWV). 



3. RESULTS 

We detect 3 mm continuum emission with a high 
signal-to-noise ratio (SNR 8-30) in all 26 sources; all 
sources remain spatially unresolved in these compact 
configuration data. Within the primary beam of ALMA 
we do not detect any other source at the sensitivity limit 
of our observations. Table [JJ lists the ALMA 3 mm con- 
tinuum positions, while the 3 mm continuum flux densi- 
ties are given in Appendix [C] together with other photo- 
metric measurements. 

Figurc[3] presents the spectra. In total, we detect 44 
line features with line integrated SNR > 5 in our survey, 
which we identify as emission lines of 12 CO, 13 CO, Ci, 
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Fig. 2. — Spectral coverage of the CO, [CI], and H2O emission 
lines as a function of redshift. The green shaded region marks 
the redshifts where two or more strong lines lines provide an un- 
ambiguous redshift, while the yellow region marks redshift range 
where only a single line is detectable. The five frequency tunings 
are shown in the left panel (see also Figure [TJ- 

H2O, and H20 + . Our spectra can be grouped into three 
categories: 

• Spectra with no line features (3 sources). 

• Spectra with a single line feature (11 sources). For 
these spectra we cannot determine the redshift un- 
ambiguously and use other spectroscopic and pho- 
tometric measurements to constrain the redshift. 

• Spectra with multiple line features (12 sources). In 
this case, a unique redshift solution can be derived 
from the ALMA 3 mm spectral scans alone. 

Table [2] summarizes the detected line features and the de- 
rived redshifts. Uncertainties for the redshifts are based 
on Gaussian fits to the line profiles. The identification of 
the ambiguous features is discussed in Section [3.21 

3.1. Additional spectroscopic observations 

For five of the sources in our sample for which we have 
detected only a single line in our 3 mm scan, we deter- 
mine the redshift using additional mm/submm or optical 
spectroscopy. We describe the observations and show the 
spectra in Appendix |A"1 

3.2. Ambiguous cases 

The most likely candidates for a single line feature 
in the 3 mm band are rcdshifted transitions of CO up 
to J=3-2 (see Figure The CO (4-3) and CO (5-4) 
lines may also appear as single lines across the band in 
cases where the Ci( 3 Pi — > 3 Po) line falls out of the cov- 
ered frequency range or may be too faint to be detected 
(the lowest flux density ratio between Ci( 3 Pi — > 3 Pn) and 
CO(4-3) or CO(5-4) that we observe in our survey is 
<0.15 (3er)). Single-line spectra cannot result from CO 
transitions of J=6-5 or higher or molecular lines that 
can a ppear at flux densities co mparable to CO (such as 
H ? O. Ivan der Werf et al.ll2010t 2011), because these lines 
would be accompanied by another line within the observ- 
ing band (see Figure [2]). The detection of FIR fine struc- 
ture lines, such as 122 ^m and 205 /im [Nil] and 158 /im 
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Source names are based on positions measured with the SPT. 
Source positions are based on the ALMA 3 mm continuum data. 

[C 11] would require extreme redshifts (z>ll) which are 
inconsistent with mm/submm continuum measurements. 

Photometric measurements allow us to discriminate 
between the possible line assignments in our single-line 
sources. The thermal dust emission of our sources is 
sampled by 3 mm ALMA, 2 & 1.4 mm SPT, and 870 /an 
LABOCA as well as 500, 350, and 250 (im Herschel- 
SPIRE observations. The photometry is given in Ap- 
pendix [D] . 

For the fitting of the therm al dust continuum w e have 
used the method described in lGreve et al.l (|2012[ ) which 
uses a greybody fit with a spectral slope of (5=2 and 
an optically thin/thick transition wavelength of 100/im, 
where the only free parameters are the d ust luminos- 
ity an d the dust temperature, Tdust- As in IGreve et al.l 
(|2012D . we exclude data points shortward of X res t = 
50/im from the fit because a single-temperature SED 
model typically cannot match both sides of the SED peak 
simultaneously due to the presence of dust at multiple 
temperatures. Both the spectral slope and transition 
wavelength affect the derived dust temperatures. For 
the present purpose, we seek only a consistent measure 
of the location of the SED peak in each source; the "tem- 
peratures" should not be interpreted as physical temper- 
atures. The dust temperature is better derived using 
the source structural information that will be available 
with lens mo dels based on high spa tial resolution ALMA 
observations (jHezaveh et al.l 120 13f) . which will help con- 
strain the dust opacity (e.g.. lWeifi et al.ll2007l ). 

Given the fundamental degeneracy between Tdust and 
rcdshift due to Wien's displacement law, it is not possi- 
ble to solve for z and Tdust simultaneously We therefore 
determine Tdust for each of the possible redshifts and 
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Fig. 3. — Spectra for target galaxies in the ALMA the 3 mm band. Spectra are shown at a resolution of 40-70 MHz (~ 100 — 250 km s 1 ) 
depending on the line width and the signal to noise ratio. 



compare these to the dust temperature distribution for 
targets with unambiguous redshifts (see Table [2]), includ- 
ing t he two SPT sources with previously known redshifts 
from lGreve et all ()2012[ ) which share the same selection 
criteria than the sample discussed here. For these sources 
we find a mean of Td U st=37.2±8.2K and no apparent 
trend with redshift (see FigureHl left). Based on the 
distribution of the temperatures in this sub-sample (19 
sources) we have calculated the probability for each of 
the four dust tcmpcrature/redshift options for the six 
sources with a single detected line and ambiguous red- 
shifts. This analysis strongly prefers a single redshift for 
one additional source (SPT0452-50, see Appendix|B|). We 
have included SPT0452-50 in the sample of sources with 
known redshifts, bringing the total to 20. 

Figure|4] shows the five remaining sources with ambigu- 
ous redshifts. We retain only dust temperatures with 
probabilities >10%, and find two plausible line identifica- 
tions/rcdshifts for four sources. In one case this threshold 
only rules out CO (2-1), leaving three plausible redshifts. 
Tablc[2] lists the possible redshifts together with the im- 
plied dust temperatures. Entries in bold face show the 



most likely redshift solution. 

In the case of the three sources without line features, 
we derive photometric redshifts based on the FIR data 
using the mean dust temperature of the objects with un- 
ambiguous redshifts. This places these three sources be- 
tween z=3.3-4.2 (sec right column of Tablc[5]). Of the two 
redshift ranges for which we cannot observe a CO line, 
the 0.36<z<l range can be excluded because the SED 
would then imply Tdust lower than the dust temperatures 
of the Milky Way and other spiral galaxies (Tdust < 15 K 
for all sources) and due to the small lensing probability 
(i j4.2[) . The galaxies may then be in the redshift desert at 
z=1.74-2.00 or at higher redshift with CO line intensities 
below our detection threshold. Our redshift survey sensi- 
tivity was intended to detect CO lines out to z~6, based 
on molecular gas estimates from the dust continuum, and 
strong detections of emission lines in 90% of the targets 
out to z — 5.7 lends credibility to the sensitivity tar- 
get. However, two of our non-detections arc among the 
1.4mm- faintest sources which leaves open the possibility 
that the line sensitivity is inadequate in these cases, al- 
beit we do detect CO lines at similar 1.4 mm flux density 
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Fig. 3. — Figurc[3] continued. 

in the survey. Yet, estimates of the CO (and C i) line in- 
tensities based on the dust continuum observations alone 
require several strong assumptions (e.g., on the gas-to- 
dust mass ratio and the molecular gas excitation). Thus 
we cannot rule out that these systems represent a class of 
galaxies with lower than expected line to continuum ra- 
tio, with the lines falling below our detection limit. If we 
place these three galaxies at z= 1.74-2. 00, we obtain low 
(Tdust~20K), but still plausible dust temperatures given 
the Tdust distribution in our sub-sample with known red- 
shifts. We note that this redshift identification is by no 
means secure, but represents the lowest plausible redshift 
range given the estimates based on the photometric data 
discussed above. 

A discussion of the 9 individual cases which have zero 
or one CO line detected with ALMA and no additional 
spectroscopic observations is presented in AppendixlCl 



4. DISCUSSION 



4.1. 



Redshift biases due to source selection criteria 

From our line identifications in Table[5J it is apparent 
that the lowest secure redshift detected in our survey is 
at 2=2.010. Only five sources are possibly at z<2 (as- 
suming that sources without a line detection fall into 
the redshift desert z=\. 74-2. 00). This is in contrast 
to the expectation from radio-identified DSFG redshift 
surveys, where typically ^50% o f all sources fall into 
the redshift range 0.5<2<2.0 (e.g.. lChapman et al.1120011 
IWardlow et~alll2011| ). 



Part of this discrepancy arises from our source se- 
lection criteria. In order to select strongly lensed, 
dusty high- redshift sources from the SPT 1.4 mm maps 
efficiently, additional criteria are used to distinguish 
the high-z population from the low-z and synchrotron- 
dominated sources that d ominate the n umbe r counts of 
Si.4mm>20mJy sources. iVieira et al.l (|2010D present a 
discussion of the classification of these populations and 
the details on how to distinguish them. Below, we pro- 
vide a summary of the selection criteria and discuss their 
impact. 

We first select sources whose mm flux is dominated by 
thermal dust emission. This step is based on the ratio 
of 1.4 to 2.0 mm flux density and is efficient at removing 
any synchrotron-dominated source from the sample, the 
majority of which are fiat-spectrum radio quasars (FS- 
RQs) and have previously been cataloged at radio wave- 
lengths. We impose a flux density cut on the sample of 
dust-dominated sources of Si. 4mm > 20 mJy based on 
the raw fluxes determined on the 1.4 mm maps. 

The second step is to use external FIR catalogs to re- 
move ('veto') low-redshift sources from the sample of 
dusty sources. Any source detected in the IR AS Faint 
Source Catalog (IRAS-FSC. fMoshir et al.llT99l at 60 or 
100/im (which implies S6o/im<200mJy and Sioo/im<l Jy 
over the entire SPT field) is omitted from our source sam- 
ple. This removes ~ 70% of the dusty sources from our 
sample. Every dusty source with a counterpart in both 
the SPT and IRAS-FSC catalogs has a published spec- 
troscopic redshift at z<0.03 and is not strongly lensed. 
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TABLE 2 
Redshifts and line identification 



Source 



T dust [K] lines 



comment 



dshifts 



secure re< 



3.0917(3) a 

4.2328(5) 

2.51480(7) 

5.699(1) 
4.2958(2) 
5.6559(4) 
4.2248(7) 
4.4771(6) 
2.0104(2) 
4.7993(5) 
3.3689(1) 
3.3988(1) 
2.1232(2) 
4.4357(6) 

4.7677(2) 
2.7799(2) 
4.5672(2) 
3.7602(3) 

2.783 
2.738 



CO(3-2 
CO (4-3 
CO(3-2 
CO(5-4 
CO (4-3 
CO(5-4 
CO(4-3 
CI(l-O) 
CO(3-2 
CI(l-O) 
CO (4-3 
CO (4-3 
CO(3-2 
CO (4-3 
CO(5-4 
CO(3-2 
CI(l-O) 
CO (4-3 

CO(7-6 
CO(7-6 



& CO(4-3) 
CI(l-O) & CO(5-4) 

& CO(6-5) 
& CO(5-4) 

CO(6-5), H 2 & H20" 1 
& CO(5-4) 
& CO(5-4) 

& CO(5-4) 
CI(l-O) & 13 CO(4-3) 
CI(l-O) & 13 CO(4-3) 

k CO(5-4) 



& CO(5-4) 
& CI(l-O) 



SPT0103-45 
SPT0113-46 
SPT0125-47 
SPT0243-49 
SPT0345-47 
SPT0346-52 
SPT0418-47 
SPT0441-46 
SPT0452-50 
SPT0459-59 
SPT0529-54 
SPT0532-50 
SPT0551-50 
SPT2103-60 
SPT2132-58 
SPT2134-50 
SPT2146-55 
SPT2147-50 

SPT0538-50 
SPT2332-53 



33.3±2.5 
31.8±3.1 
40.7±4.2 
30.1±4.9 
52.1±7.8 
52.9±5.3 
52.9±7.5 
39.3±3.9 
20.9±1.8 
36.0±3.7 
31.9±2.4 
35.1±3.0 
26.3±2.0 
38.6±3.5 
37.8±4.5 
40.5±4.6 
38.7±5.1 
41.8±4.1 

31.2±7.1 
32.9±3.6 



3.9592(5) 

2.7174(6) 
3.5956(3) 

2.4474(3) 
3.6854(2) 

4.8565(2) 

2.5142(1) 
2.2335(2) 

1.1557(1) 
3.1286(5) 

2.0966(4) 



CO(4-3 
CO(3-2 
CO (4-3 
CO(3-2 
CO(4-3 
CO(5-4 
CO(3-2 
CO(3-2 
CO(2-l 
CO(4-3 
CO(3-2 



CO(l-O) from the ATCA 



CI(l-O) feature low SNR, [Cll] confirmation with APEX 
alternative redshifts excluded due to lack of higher J CO lines 



VLT Civ 1550 A detection 



[CII] from APEX 

CO(7-6) & CO(8-7) detections from Z-Spec 



the SMA 



CO(8-7), Si iv 1400 A ZSpcc/VLT from Greve et al. (2012) ; no ALMA data 
Lyq, Civ 1549 A ZSpcc/VLT from Greve et al. (2012); no ALMA data 

ambiguous redshifts 



CI(l-O) 
CI(l-O) 



CI(l-O) feature low SNR 
alternative ID if CI(l-O) is not real 
CI(l-O) feature low SNR 
alternative ID if CI(l-O) is not real 

similarly likely ID 



SPT0125-50 
SPT0300-46 
SPT0459-58 

SPT0512-59 
SPT0550-53 



43.3±5.2 
29.5±3.2 
38.6±3.6 
26.7±2.2 
32.0±4.5 
40.8±6.0 
22.4±2.9 
33.2±3.0 
20.4±1.6 
30.6±4.6 
21.6±2.9 



no CO line detections 



SPT0128-51 - - no lines 2=1.74-2.00 ?; z photo = 3.8 ± 0.5 for T dust =37.2 K 

SPT0319-47 - no lines z=l. 74-2.00 ? ; z photo = 4.2 ± 0.2 for T dust =37.2K 

SPT0457-49 - no lines z=1.74-2.00 ? ; z photo = 3.3 ± 0.2 for T dust =37.2K 



Note. — In case of ambiguous redshifts, preferred solutions are shown in bold. 
a The number in brackets is the redshift uncertainty in the last decimal derived from Gaussian fits to the line profiles. 



The third step is to use external radio catalogs to 
remove low-redshift and radio-loud sources from the 
sample of dusty sources. Any source detected in 
the 843 MHz Sydney University Molonglo Sky Survey 
(SUMSS. iBock et al.llilfflSl) (with a -6mJy 5cr flux den- 
sity threshold over the entire SPT field) is omitted from 
our source sample. The SUMSS veto removes an ad- 
ditional ~ 15% of the dusty sources which passed the 
IRAS veto. This step is intended to ensure that no FS- 
RQs were allowed into the sample. The mean radio flux 
density reported in the SUMSS catalog for these sources 
is (S , 843MHz)=52mJy, well above the catalog threshold. 

The effect of these selections on the rcdshift distribu- 
tion of the 1.4 mm sources targeted in this study depends 
on the intrinsic radio-IR SEDs of the DSFGs. Figurc[5] 
shows the redshift limits beyond which different radio- 
IR SEDs pass our source veto criteria. We show here 
well-studied examples of quiescent and star-forming local 
galaxies, as well as an example for a high- redshift, radio- 
loud active galactic nucleus (AGN) host galaxy. The 
figure demonstrates that galaxies which follow the local 
radio-FIR correlation and have relatively cold dust tem- 



peratures (Tdust^$30 K, e.g., M51) would pass our source 
selection criteria at relatively low redshift (z>0.5). 

Sources with Arp 220-like SEDs would pass our selec- 
tion criteria at higher redshifts (z>1.4). Other local and 
high-z IR luminous so urces, including M82, SMM J2135- 
0102 ('The Eyelas h' - iSwinbank etHfll 121)10) . and HR10 
(|Stern et al.l 120061 not shown), are allowed at redshifts 
similar to Arp 220. Sources with FIR SEDs dominated 
by hotter dust (due to AGN h eating, as in H1413 +117, 
also known as 'The Cloverleaf': lBenford et al1 ll999) than 
is typical for star-forming systems would be found in 
IRAS and excluded from the sample out to z ~ 3. 

The SUMSS veto may exclude a few source classes 
from our sample. Figure [5] shows that systems with 
much higher radio power than implied by the radio-IR 
correlation, such as lensed radio-loud AGN with signif- 
icant dust emission (e.g., the Clovcrlcaf), are excluded 
from our sample over a large redshift range. This veto 
may also exclude lensed DSFGs at z<1.5 (coincidentally 
close to the IRAS redshift veto limit), where the radio- 
FIR correlation predicts the radio emission will exceed 
the SUMSS limit. Finally, DSFGs lensed by foreground 
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Fig. 4. — Le/t: Dust temperatures for the sources in our sample with unambiguous redshifts. Center: Combined histogram of dust 
temperatures derived from the posterior likelihood distributions for the sources with unambiguous redshifts. Ovcrplotted are the dust 
temperatures determined for each redshift option for those sources with uncertain redshift; horizontal spacing is arbitrary. The solid 
and dashed lines show the median and 95% confidence interval dust temperatures for those sources with unambiguous redshifts. Right: 
Probability for the single line detected in our ALMA spectrum to be identified as one of the four possible CO transitions for the five sources 
with ambiguous redshifts. The probabilities were calculated by comparing the dust temperature associated with each line identification to 
the dust temperature distribution of our sources with known redshifts. The horizontal dashed line shows a probability of 10%, the cut off 
above which we consider the line identification to be plausible. 

galaxies with radio-active AGN and residual FSRQs will 
be excluded in a redshift-unbiased way by this veto. 
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Fig. 5. — Top: Redshift bias due to our IRAS 60 and 100/im, 
and 843 MHz radio flux vetos. The bars show the redshift range 
for which specific radio-to-IR SEDs arc excluded from our sample. 
The color coding of the bars corresponds to galaxies shown in the 
bottom part of the figure. Bottom: Radio to optical SEDs of M51 
(the Whirlpool Galaxy), Arp220 (the nearest ultraluminous in- 
frared galaxy) and H1413+117 (the Cloverleaf QSO). These galax- 
ies represent a range of possible SED types and are normalized 
to Si. 4 mm = 28mjy (the mean 1.4mm flux density of our sam- 
ple). The dashed horizontal line shows our selection wavelength 
of 1.4mm. The arrows show the 843 MHz, lOOfim and 60/im up- 
per limits used for our source selection. The SEDs are shown for 
the lowest redshift (value indicated in the figure) for which each 
source matches our selection criteria, except for H1413+117 which 
is shown at z = 3.0. 




Fig. 6. — Probability of strong gravitational lensing as a func- 
tion of redshift fo r different source m agnifications (fi) calculated 
from the models of Hczavch & Holder (2011) . The model assumes 
no size evolution for the underlying DSFG population. The fig- 
ure demonstrates the strong decline of the lensing probability for 
z<1.5, independent of the magnification. 



4.2. Redshift biases due to gravitational lensing 

The high 1.4 mm flux density cut of our target selection 
implies that even the most infrared-luminous galaxies are 
too faint to be included in the SPT dusty-source sample 
at z > 0.5 without assistance from gravitational lensing 
(Lib. > 3 • 10 13 L Q for a Arp220 like SED). This expecta- 
tion is confirmed by our ALMA high-angular resolution 
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imaging that resolves our sources into arcs or Einstein 
rings - hallmarks of gravitational lensing (|Vieira et al.l 
[2013h . The rcdshift-dcpcndent probability of strong 
gravitational lensing therefore has important effects on 
our redshift distribution. In Figure[6j we show the differ- 
ential probability of stron g lensing versus red s hift, c alcu- 
lated from the mode ls of iHezaveh fc Holderl (|2011[ ) and 
iHezaveh et al.l (|2012f) , which use gravitational lensing by 
a realistic population of intervening halos to match the 
observed number counts of bright DSFGs. The strong 
evolution in the lensing probability (the fractional vol- 
ume at each redshift subject to high magnification), a 
factor of 20 between z~2 and z ~ 0.5, demonstrates that 
the requirement that we find lensed sources strongly sup- 
presses sources at z<1.5. For z>2 the lensing probability 
varies much more slowly, implying weaker effects on the 
lensed source counts. 

At higher rcdshifts, other lensing effects can more 
significantly alter the normalized redshift distribution, 
dn/dz, especially changes in source sizes. To evaluate 
such effects, we compare an assumed intrinsic redshift 
distribution to the model distribution of strongly lensed 
sources (Si.4 mm >f 5 mJy, consistent with the deboosted 
1.4 mm fl ux densities of our sou rces, see Table |H). As dis- 
cussed in IHezaveh et al.l (|2012l ).the selection of a sample 
of millimeter-bright DSFGs, lensed by intervening galax- 
ies, will preferentially identify those with more compact 
emission regions. This implies that the observed redshift 
distribution could be biased if DSFGs undergo a size evo- 
lution with redshift. 

Observationally, it is well established that high-rcdshift 
DSFGs are significantly larger than local ULIRGs. In 
the high-rcdshift (z>2) sources the star-forming regions 
extend over ^5 kpc diameter, while lower-redshift (z<l) 
ULI RGs typically f orm stars in kpc-sized r egions (see, 
e.g., iTacconi et all 120061 lEngel et al.l 120101 and refer- 
ences therein). Whether DSFGs undergo a size evolu- 
tion in the redshift range z=1.5-6, the relevant redshift 
range for our study, is, however, largely unknown due 
to the small number of high-redshift objects for which 
spatially resolved observations of the submm emission 
region exist and the large diversity of morphologies. Ev- 
idence for extended molecular gas reservoirs (>10kpc 
diameter) has b een found in some DSFGs out to red- 
shift z«4 (e.g . iGenzel et al.l l200l llvison et al.l 120101; 
Younger etalT 12010] iCarilli et al.l 120 lit llvison et al.l 



201 It Riecher s et al while the molecular gas dis- 

tribution in IR luminous AGN host galaxies, which have 
been measured out to redshift ,2=6.4 , are typically more 
comp act (^2-3 kpc diameter, e.g., IWalter et al.l |2004[ 
120091) . These differences, however, mainly reflect the di- 
versity of submm-detected objects and possibly an evo- 
lutionary link betwee n DSFGs and AGN host galaxies 
(jRiechers et al.ll2011al ) rather than an overall size evolu- 
tion of submm-selected high-z galaxies. 

In Figure [JJ we compare different size-evolution sce- 
narios, where the intrinsic distribution was prescribed 
to be consistent with the observed redshift distribution 
from radio-identified DSFGs inc luding recent spectro- 
scopic data f rom t h e literature (IChapman et al.l 120051: 
Capak etafl 120081: iCoppin et al.l 120091 iDaddi etal l 
1.1 120101: If" 



2009allbl: iRiechers et al.1 120101: IBanerii et al.1 120111: 
Walter et al.l 120121 ^ The figure demonstrates that the 



effect of gravitational lensing on the observed redshift 



distribution is relatively small when there is no size 
evolution or increasing source sizes with redshift. For 
example, in the redshift range z=2-4 the difference 
between dn/dz derived from the unlensed and lensed 
sources is smaller than ~ 20% in both cases. In the case 
of no size evolution the observed redshifted distribution 
is displaced by Az^0.3 towards higher redshifts com- 
pared to the unlensed case. Given the steep increase 
of dn/dz betwe e n z= l -2 of the redshift distribution 
([Chapman et al.l 120051 : IBanerii et all 120111 ). this shift 
causes an underestimate of the source counts in this 
redshift interval by roughly a factor of two which may 
explain the low number of z<2 objects detected in our 
survey. For decreasing source sizes with redsh i ft (as 
suggested by optical observations, iFathi et al.1 120121) 
the difference between the observed and intrinsic red- 
shift distribution can become significant also for z>3, 
with the counts of the high-redshift galaxies increased 
compared to the intrinsic distribution. 

A compilation of the effective source radii for z=l-6 
derived from an analysis of the dust SEDs of unlensed 
submm detected DSFGs and quasi- stellar object (QSO ) 
host galaxies has been published in iGreve et al.1 (|2012D . 
Their Figure 5 shows the submm source radii as a func- 
tion of redshift. The size of the highest redshift sources 
(z=5-6) in this diagram tend to fall below the average 
size of z=l-3 objects, but as mentioned above, these 
high-redshift sources are all QSO host galaxies and as 
such cannot be taken as evidence for a size evolution of 
the whole DSFG population. Th e sa mple of source radii 
in the literature (ITacconi et al.l 120081: lEngel et alJl201Ct 
iRuiopakarn et aLl^OUl T which were directly measured 
from high-resolution imaging, show no clear evidence for 
size evolution above z>0.4. In the absence of conclu- 
sive observational constraints, it is difficult to quantify 
the redshift bias due to gravitational lensing. We note, 
however, that making our observed redshift distribution 
consistent with a n intr insic distribution like the one from 
IChapman et al.l (|2005l ) would require an extreme growth 
of DSFGs between z=6 and z=2 (r=0.2kpc to 2.5 kpc 
in 2.3 Gyr, see Figure[7J. Likewise, a modest evolution 
(r=1.5kpc at z=6 to 2.5 kpc at z—2, using the QSO size 
measurements as lower limits to the size of DSFGs at 
z = 6, see above) results in a steeper redshift distribu- 
tion than that implied by our most likely redshifts. Both 
suggest that gravitational lensing is unlikely to be the 
dominant source for the differences in dn/dz between 
the present sample and the radio-identified samples. 

4.3. The redshift distribution 

Even with the conservative choice of taking all ambigu- 
ous sources to be at their lowest redshift option (see Table 
[2]), at least 50% of the SPT sample is at z>3. Only five 
sources are possibly at z<2 (assuming that sources with- 
out a line detection fall into the redshift desert z=1.74- 
2.00), consistent with the expectations for a sample of 
strongly lensed objects. Our sample mean redshift is 
z=3.5. This redshift distribution is in contrast to that of 
radio identified DSFGs which have a significantly lower 
mean redshift of z=2.3 and for which on ly 10-15% of the 
popu l ation is expected to b e at z>3 (e.g. . IChapman et al.l 
l2^CMrWardlow et al.ll20ll . 

A potential difference between our redshift distribution 
and the 850-/mi-selected samples in the literature arises 
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TABLE 3 

Measured redshift distribution for SPT 
sources 



Fig. 7. — Comparison between an assumed intrinsic redshift dis- 
tribution (dn/dz, solid black) consistent with spectroscopic obser- 
vations (see text for references) and distribution s modified by grav- 
itation al lcnsing (using the models described injHczaveh fc Holder! 
Q2011f l) under different size evolution scenarios. The green dashed 
line shows the bias to the redshift distribution due to gravitational 
lensing assuming no size evolution versus redshift. The blue dot- 
dashed line show the bias to the redshift distribution due to grav- 
itational lcnsing if the size of DSFGs increases with redshift, from 
r=l kpc at z=2 to 3 kpc at 2=5. The red dotted line shows the bias 
of the redshift distribution due to gravitational lensing if the size 
of DSFGs decreases moderately with re dshift (r=2.5kpc at 2=2 as 
measured for DSFGs (Engcl ct al. 201( j), to r=1.5kpc at 2 =6 using 
the measured submm QSO host sizes TWalter et aUl2009h as lower 
bound to the size of DSFGs). The maroon dotted line exemplifies 
the extreme size evol ution which would be r equired to bring the 
redshift distribution of Chapman ct al. (2005) into agreement with 
our observations (r=2.5kpc at 2=2 to r=0.2kpc at 2=6). 

from the interaction of the SED of the typical DSFG 
and the selection wavelength. This has be e n disc ussed 
in several papers, including iGreve et al.l (|2008[ ) and 
ISmolcic et al.l (|2012t) . It has been argued that 850 /im 
selection results in lower redshift samples than 1.4 mm 
selection because the negative if-correction ceases once 
the SED peak is redshifted into the detection band, which 
occurs at lower redshift for shorter wavelength observa- 
tions. Because our sources have been selected at 1.4 mm 
(SPT) and also observed at 870 urn (LABOCA), we can 
examine the effect that 850 jim selection would have on 
our sample. The flux ratio as a function of redshift is 
shown in FigureEl it reveals a modest decrease of the 
870/im/1.4mm flux ratio for increasing redshift. Our 
observations therefore support the notion that 850 fim 
selection will preferentially remove sources at the high- 
est redshifts. We caution, however, that this effect will 
operate only on the fainter population of high-rcdshift 
sources, those near to the detection limit where the 
850 /im may fall below the detection threshold while the 
1.4 mm signal remains detectable. 

Some studies of submm selected galaxies from blank 
field surveys presented evidence for a correlation between 
observed submm flux density and the source redshift (e.g . 
Ivison et al l 120021 : iPope et all 120051 : llvison et al.1 120071 : 
Biggs et al.ll2011 ). If confirmed, this could imply a pos- 
sible bias towards higher redshift for our study if the in- 
trinsic IR luminosity of our sample is on average higher 
than that of unlcnsed mm/submm selected samples. So 
far, lens models based on spatially resolved images of 
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the 870/ Ltm continuum are o nly available for four SPT 
sources (jHezaveh et al.ll2~013l ) . These have magnifications 
of /i = 5—21 with a mean of /2=14. The gravitational flux 
a mplification o f the S PT sources has also been discussed 
in IGreve et"aD (|2012t ). They derive p, = 11 - 22 based 
on an analysis of the FIR properties of 11 SPT sources 
compared to unlensed samples, in reasonable agreement 
with the lens models. Adopting an average magnifica- 
tion of fi — 15 for the sources studied here, our sample 
is expected to cover intrinsic flux densities of Si.4 mm = 
1.0 — 3.0 mjy and Ssro^m = 1-7 — 9.5 mJy with means 
of Si. 4 mm = 1.8 mJy and Ss70pm = 5.4 mJy. These in- 
tensities ranges arc well match with unlensed source flux 
densi ties observed in mm/submm blank fi e lds surveys 
'e.g iBorvs eFail [2001 jCoppin et all [200l Wove et al.l 



[2Tj0aTAl!ist"ermann et alj|2009t iWeifl' et al. [ |2009bD which 
implies that our sample should be representative for the 
submm selected galaxy population at z > 1.5. We fur- 
ther note that the claimed correlation between observed 
submm flux density and source redshift has recently been 
questioned (|Wardlow et alJl20lH Ikarim et al.ll2012r i. 

An additional difference between this sample and ear- 
lier spectroscopic measurements of the DSFG redshift 
distribution is the radio selection. As noted above, pre- 
vious DSFG redshift searches have primarily relied upon 
radio counterpart identification to provide optical spec- 
troscopy targets and therefore have a radio detection re- 
quirement. Here we have excluded sources with bright 
radio counterparts, which might be expected to oppo- 
sitely bias the sample. However, a comparison of the 
submm-radio flux d ensity ratio distr i bution for the radio- 
identified sample of lChapman et ail (|2005l ) and the simi- 
lar ratio (corrected for differences in observing frequency) 
constructed from SUMSS and SPT measurements for our 
SUMSS-vctocd sources shows that these objects emit a 
much larger fraction of their en ergy in the radio than 
even the most extreme sources in I Chapman et all ()2005l ) 
(see their Figure 7). Likewise, sources that pass our 
SUMSS radio- veto are not biased towards larger submm- 
radio flux density ratios than radio selected samples from 
the literature due to the shallowness of the SUMSS sur- 
vey. Therefore this veto should not preferentially exclude 
low-redshift DSFGs, though optical spectroscopic mea- 
surements of the excluded sources will be useful in de- 
termining which source classes and which redshifts dom- 
inate the excluded objects. 

The determination of the shape of our redshift dis- 
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Fig. 8. — Observed 870 fim to 1.4 mm flux density ratio as a func- 
tion of rcdshift for our sample of 20 sources with unambiguous spec- 
troscopic rcdshifts. The grey points show the individual measure- 
ments and their error bars taking absolute calibration uncertainties 
into account. The black crosses show the mean flux density ratio in 
rcdshift bins of Az=l centered at the weighted mean z. The dashed 
line is a linear fit to the data (Sg70fim/ Si. 4mm = 4.18 — 0.34z for 
2 = 2 — 6). The dotted line shows the expectation for a Arp220 
like dust SED. 

tribution is currently hampered by the eight ambiguous 
rcdshifts. In Figurc|9] (left) wc compare two redshift dis- 
tributions, one using the lowest redshift option for all 
sources, and the other assuming the most likely redshift. 
In the first case, our rcdshift distribution shows some 
evidence for a peak at z « 3, consistent with the find- 
ings of radio identified DSFGs, and then decreases out 
to z^6. The decrease, however, is much shallower than 
suggested from radio identified DSFGs. In the latter case 
our redshift distribution rises up to z ~ 4 and falls off at 
higher rcdshift. Within the errors both distributions are 
consistent with a fiat redshift distribution between z~2- 
4. Note that to these distributions we hav e added two 
additi onal strongly lenscd SPT sources from lGreve et al.l 

(pil . 

We adopt the redshift distribution informed by our 
dust temperatures and other data ("SPT best" in Fig- 
ure [9]) for the discussion which follows, and report the 
values for dn/dz in Table 

FigureO (center) highlights the large difference be- 
tween our results and previous redshift surveys. 
Compared to previous surveys with spectroscopic 
rcdshifts that rely on radio counterpart identi- 
fication (IChapman et all 120051 iCapak et al.l 120081: 



Coppin et al.l 120091 iDaddi et al.ll2009aObt iRiechers et al.1 
20101 : iBanerji et alJl20TH iWalter et aT]|2012f > we find a 
far greater fraction of high-redshift sources. As discussed 
earlier, gravitational lensing may explain part of this dis- 
crepancy if DSFGs are smaller at high redshifts, though 
extreme evolution is required to explain the full differ- 
ence. Recent work based on CO(l-0)-derived redshifts 
for a DSFG sample selected from the H-ATLAS survey 
([Harris et al.l 120121 not shown here) implies a r edshif t 
distribution in agreement with IChapman et al.l (|2005l ). 
These sources, however, were selected to peak in the 
SPIRE 350/im channel to match the 2.1<z< 3.5 red- 
shift coverage of the instrument used to measure rcdshifts 



(jHarris et al.l [2012). Despite this selection, > 50% of 
their targeted sources remained undetected in CO, which 
may imply that there are a significant number of sources 
a t redshifts larg e r than z=3.5 in this sample as well. 

iSmolcic et al.l f|20 12f ) also find an increased fraction of 
DSFGs at z>3 through a combination of spectroscopic 
and photometric redshifts for a mixed sample of 1.1mm 
and 870/mi selected sources in the COSMOS field. They 
note that 50-70% of their z>3 DSFGs have no radio 
counterpart down to ~10 /iJy at 1.4 GHz, which supports 
the prediction that including radio counterpart identifi- 
cation in the process of surveying DSFG redshifts will 
suppress highcr-z sources, as expected from SED tem- 
plates. The similarity in t he redshift distributi on of un- 
lcnscd sources compiled by ISmolcic et al.l (I2012D . derived 
primarily from photometric redshifts, and our own (Fig- 
ure [9j center) may be evidence that gravitational lensing 
is not strongly affecting the underlying rcdshift distribu- 
tion. However, greater numbers of molecular-line-derived 
rcdshifts for both populations will likely be required to 
settle this issue. 

In the case of no size evolution in DSFGs, our study 
suggests that previous spectroscopic DSFG redshift sur- 
veys, which are almost exclusively based on radio iden- 
tified sources, have missed >50% of the DSFG popula- 
tion as it resides at redshifts z>3 and the putative high- 
redshift tail of DSFGs may in fact turn out to be a much 
broader, flat-topped rcdshift distribution which could ex- 
tend to z > 4. 

4.4. Comparison to models 

Redshift distributions (dn/dz) and number counts are 
the main observational constraints to galaxy formation 
models. Matching available data for DSF Gs with these 
mode ls has been particularly difficult (e.g.. iBaugh et all 
120051 ). requiring some ad hoc changes such as top-heavy 
initial mass functions. As argued above, our dn/dz - 
although currently based on only 28 sources — appears 
significantly different from the current ly largest sam- 
)le of spectroscopic DSFG rcdshifts by IChapman et al.l 
2005h . With direct mm identifications, a 71% spectro- 
scopic completeness, and likely redshifts for an additional 
18%, our SPT DSFG dn/dz represents an important new 
observational constraint to these models. 

We compare our measured dn/dz with four recent 
models in Figure [9] (right), removing sources at z<1.5 
from the models to mimic the strong lensing selection 
selection described in Section l4~3l We discuss the individ- 
ual models below and give the % 2 for each model for the 
five redshift bins. Despite the relatively small number 
of redshifts, our new SPT dn/dz already discriminates 
b etween galaxy formation models. 

iBethermin et al.l (|2012f ) present an empirical model 
starting from the observed FIR number counts split into 
"main-sequence" and starburst mode star-forming galax- 
ies. Their model includes the effects of magnification 
by strong lensing, so it can directly predict the dn/dz 
for the SPT sample. For the comparison with our data 
we use the predicted redshift distribution for sources 
with Si.4 mm >15 mJy, consistent with our source selec- 
tion. This model matches our rcdshift histogram very 
well, with a comparison to the five rcdshift bins giving a 
X 2 of 1.9 across five reds hift bins. 

The lLacev et al.l (|2010f ) model is a semi-analytic model 
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Fig. 9. — Left: Redshift distribution of strongly lensed DSFGs derived from our molecular line survey. The red histogram shows the 
2-distributions for the SPT sources using the most likely redshift redshift identification for the sources with ambiguous redshifts; the 
blue histogram shows the same for the lowest redshift identification o f these five sour c es (se e TablpO. Middle: Redsh ift distrib u tion o f 
radio-identified DSF G s with spectros c opic r edshifts at z> 1.5 (g r een, IChap man ct al] i2005ft ; [Capak ct al. (200|); ICoppin et all ((2009); 
IDaddi et all 12009a b): Rlc chers et"al1 MOffl ); IBanerii etUI (l2TmT ) ; IWalter et al.l 12013) ). mm-identified DSFGs with photometric redshifts 
at z>1.5 from (black. ISmolcic et al.l 112012 1 )). compared to the mo st likely SPT dis t ributi on (Red). Right: Red s hift d istributions from the 
models of IHavward et al.l <2012F 7BlackT lBetherrniii et al.l 12012) (Blue), IBensonl 120121) (Green) iLacev et al.l 120101) (Purple) for z>1.5, 
compared to the most likely SPT distribution (Red). The arrows in all panels show the mean redshift of each distribution. In all panels 
the histograms are calculated over the same redshift bins but are plotted with slight shifts in z for clarity. 

methods typically used to derive DSFG redshifts. 

In total we detect 44 spectral features in our survey 
which we identify as redshiftcd emission lines of 12 CO, 
13 CO, C i, H 2 0, and H 2 + . We find one or more lines in 
23 sources, yielding an unprecedented ~90% success rate 
of this survey. In 12 sources we detect multiple lines. In 
11 sources we robustly detect a single line, and in one 
of those cases we can use that single line to obtain an 
unambiguous redshift. For an additional five galaxies, in 
which we detect a single line with ALMA, we can deter- 
mine the redshift using additional spectral and optical 
data yielding 18 unambiguous redshifts. For five sources 
with a single line detection we have used our excellent 
mm/submm photometric coverage (3 mm to 250 fim) to 
narrow the line identification and make a probabilistic 
estimate for the redshift based on the FIR dust tem- 
perature derived from extensive broad band photometric 
data. In three sources we do not detect a line feature, 
either because the lines are too weak, or because they are 
in the redshift desert z~ 1.74-2. 00. Adding in two previ- 
ously reported SPT sou rces with spectroscopic redshifts 
fromjG reve "et~aT1l20rl . we derive a redshift distribution 
from 28 SPT sources. 

We analyze the redshift biases inherent to our source 
selection and to gravitational galaxy-galaxy lensing. Our 
selection of bright 1.4 mm sources imposes a requirement 
that they be gravitationally lensed, effectively suppress- 
ing sources at z < 1.5 due to the low probability of being 
lensed at these redshifts. Beyond z ~ 2, gravitational 
lensing does not significantly bias the redshift distribu- 
tion unless DSFGs undergo a systematic size evolution 
between z=2-6 with decreasing source sizes for higher 
redshifts. An analysis of the black body radii of un- 
lensed DSFGs from the literature does not support the 
existence of such an evolution, but it also cannot be ex- 
cluded conclusively at this point. 

Our sample mean redshift is z=3.5. This finding is in 
contrast to the redshift distribution of radio identified 
DSFGs which have a significantly lower mean redshift 
of z=2.3, and for which only 10 -15% of the populatio n 
is expected to be at z>3 (e.g., lChapma"n 
The redshift distribution of our sample appears almost 
flat between 2; =2-4. Our study suggests that previous 



identical to that presented in iBaugh et al.l (|2005l) . The 
model employs a top-heavy stellar initial mass function, 
which results in more luminosity and more dust produced 
per unit star formation rate, to better match the bright 
end of 850 /im galaxy counts. This model does not in- 
clude the effects of strong lensing, and DSFG counts are 
based on a selection in Si.4 mm with > 1 mJy (C. Lacey, 
private communication). The \ 2 between this model and 
our m easurement acro ss the five redshift bins is 10.7. 

The IBensonl (|2012f ) model is a semi-analytic mode l 
that also exp ands upon the work of IBaugh et al.l (|2005h . 
Whereas the ILacev et al.l (|2010l ) model required a top - 
heavy stellar initial mass function, the IBensonl (|2012f ) 
model merely has enhanced dust production in star- 
bursts. This model does not include the effects of strong 
lensing, and DSFG counts are based on a selection in 
Ssso/xm (> 5 mJy) The predicted dn/d z distribution 
comes close to the IChapman et al. l (120051) distri bution, 
but clearly fails to fit the SPT or ISmolcic et al.1 (|2012f ) 
measurements. Part of this difference may be due to the 
850 /an instead of 1.4 mm source selection, and a pos- 
sible lensing bias. The x 2 between this model and our 
measurement across the five redshift bins is 39.8. Our 
measurements a re clearly at odds with this model. 

The model by IHavward et ail (|2012l ) combines a semi- 
empirical model with 3D hydrodynamical simulations 
and a 3D dust radiative transfer. Strong lensing is not in- 
cluded in the modeling and the model predicted dn/dz is 
determined using sources with 5i mm >lmJy, consistent 
with the expected intrinsic flux densities of our sample. 
The distribution of the DSFGs in this model is close to 
the observed SPT dn/dz, with a x 2 of 2.8 between data 
and model. 

5. SUMMARY AND CONCLUSION 

We have used ALMA to measure or constrain the red- 
shifts of 26 strongly lensed DSFGs detected in the SPT- 
SZ survey data. The redshifts were derived using molecu- 
lar emission lines detected in frequency scans in the 3 mm 
transmission window covering 84.2 to 114.9 GHz. As the 
molecular emission lines can unambiguously be associ- 
ated with the thermal dust continuum emission at our 
selection wavelength of 1.4 mm, this technique does not 
require any multi-wavelength identification unlike other 
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spectroscopic redshift surveys of DSFGs based on radio 
identified sources are likely biased towards lower redshift 
and have missed a large fraction (> 50%) of the DSFG 
population at rcdshifts z>3. 

With a 90% detection rate, our ALMA+SPT CO red- 
shift survey is the most complete DSFG survey to date. 
It demonstrates the power of ALMA, with its broad-band 
receivers and large collecting area, to provide the critical 
galaxy redshift information needed to measure the cos- 
mic history of obscured star formation, particularly at 
the highest redshifts where other techniques falter. The 
magnification of the S PT sources by inter vening mass 
(factors of ~10 or more, Hc zaveh et al.ll20~13l ) has allowed 
us to obtain these results in the early science phase of 
ALMA, with only 16, of the eventual array of 54, 12- 
meter antennas. With the full array, such studies will 
be possible on unlensed sources, highlighting the enor- 
mous scientific impact ALMA will have in the coming 
decades. With spectroscopic redshifts for a large num- 
ber of DSFGs, it is now possible to study the conditions 
of the interstellar medium at high redshift in great detail 
through spatially resolved spectroscopy of FIR molecular 
and atomic lines. The SPT sources presented here repre- 
sent less than 25% of the entire sample of high-redshift, 
strongly lensed DSFGs. Obtaining redshifts for the re- 
maining sources will enable us to definitively constrain 
the rcdshift evolution of DSFGs. 
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APPENDIX 

SUPPLEMENTARY REDSHIFT INFORMATION 

In this appendix, we show the supplementary observations that resolve redshift ambiguities in our ALMA observa- 
tions: 

SPT0125-47: The identification of the 98 GHz line as CO(3-2) is confirmed with a CO(l-O) detection using the 
Australia Telescope Compact Array (Figure fTUj) . 

SPT0441-46'- The identification of the 105 GHz line as CO(5-4) is confirmed with a [CII] 158^m detection with 
the First Light APEX Submillimetre Heterodyne receiver (FLASH) on APEX (Figure [Til) ■ The low S/N [CI] (1-0) 
detection with ALMA further strengthens this redshift identification. 

SPT0551-50: A strong emission line is visible at ^4800 A using the VLT FOcal Reducer and Spectrograph (FORS2; 
lAppenzeller et aL1ll998D . which is consistent with the 3 mm CO (3-2) line if we ascribe it to CIV 1550 A. See Figure [T2l 

SPT2134-50 : The CO(7-6) and C O(8-7) lines are detected in a 190-310 GHz spectrum (Figure H3J) obtained with 
Z-Spcc/APEX (Bradford et al. 2004), and subsequently confirmed through Submillimeter Array (SMA) observations 
of CO(7-6) and [CI] (2-1) (See Figure[H]). The ALMA data, released later, agree with this identification, with ALMA 
detecting the CO(3-2) line at 91.5 GHz. 

SPT2132-58: The identification of the 100 GHz line as CO(5-4) is confirmed with a [CII] 158/mi detection with the 
First Light APEX Submillimetre Heterodyne receiver (FLASH) on APEX (Figure fT5|). 

SUPPLEMENTARY INFORMATION FOR SOURCES WITH A NO OR SINGLE LINE DETECTIONS 

Below, we discuss the 9 individual cases which have zero or one CO line detected with ALMA and no additional 
spectroscopic observations. 
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Fig. 10. — Australia Telescope Compact Array spectrum of SPT 0125-47 showing the CO(1-0) line confirming the single ALMA line as 
CO(3-2) at 2=2.5148. 
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Fig. 11. — APEX/FLASH spectrum of SPT 0441-46 showing the [CII] A158£tm line (filled histogram) confirming the single ALMA line 
as CO(5-4) (red line, scaled to allow for a comparison between the line profiles) at 2=4.4771. 
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Fig. 12.— VLT/FORS2 spectrum of SPT 0551-50 showing the CIV A1549A line confirming the single ALMA line as CO(3-2) at 2=2.123. 
Thin red dashed lines indicate the wavelengths of expected spectroscopic features, while thick green dotted lines mark areas dominated by 
skylines. 



SPT0125-50: In this galaxy we detect a second tentative line feature at 99.20 GHz which is consistent with the 
expected frequency for Ci( 3 Pi — ^ 3 Po) if the 93.03 GHz line is CO(4-3). This is our preferred identification, giving 
z=3.959. In case the weak 99.20 GHz feature is not real, CO(5-4) as identification for the bright line can be excluded 
as CO(6-5) should have been detected too. For CO(2-l) at z=1.343, the implied dust temperature would be 17 K, 
lower than any we observe. An additional plausible identification is CO(3-2) at 2=2.717 (Td U st=30K). 
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Fig. 13.— APEX/Z-spec spectrum of SPT 2134-50 showing 2-3ct detections of the CO(7-6) and CO(8-7) lines confirming the single 
ALMA line as CO(3-2) at 2=2.779. Dashed lines mark the expected frequencies of CO and H2O features. The combined significance of 
the lines detections is 5. 60". 
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Fig. 14. — SMA spectrum (filled histogram) of SPT 2134-50 showing CO(7-6) and evidence for [CI] (2-1) confirming the single ALMA 
line as CO(3-2) (red line, scaled to allow for a comparison between the line profiles) at 2=2.779. 
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Fig. 15.— APEX/FLASH spectrum of SPT 2132-58 showing the [CII] A158£tm line (filled histogram) confirming the single ALMA line 
as CO(5-4) (red line, scaled to allow for a comparison between the line profiles) at 2=4.7677. 
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SPT0128-51: No line is detected in this spectrum. If it is in the 2=1. 74-2. 00 rcdshift desert, the dust temperature 
is a low Tdust ~ 19 K. Alternatively, at higher rcdshift the line-to-continuum ratio should be smaller and could go 
undetected. If SPT0128-51 has the same Tdust as the median temperature of the unambiguously identified population, 
37 K, its corresponding photometric redshift would be z — 4.3. 

SPT0300-46: This source is similar to SPT0125-50 and has a clear CO detection at 100.30 GHz and a tentative 
Ci( 3 Pi -+ 3 P ) line at 107.08 GHz which implies CO (4-3) at 2=3.594. If the latter feature is not real, CO(3-2) at 
2=2.446 and Td U st=27K is an alternative interpretation. CO(2-l) at 2=1.298 would imply Td U st=17K, which we 
consider unlikely. CO(5-4) can be ruled out as CO(6-5) would have also been detected. 

SPT0319-47: No line is detected in this spectrum. The dust temperature would be « 20 K if the source is in the 
2=1.74-2.00 redshift desert. As with SPT0128-51, a higher redshift with weak lines cannot be ruled out. Matching 
this source to the median temperature of the known sample yields a photometric redshift of z = 4.0. 

SPT0452-50 : There is a clear line detection at the very edge of the band (114.87GHz). CO(4-3) and CO(5-4) can 
be excluded as a second CO line would be detected in the band. CO(2-l) at 2=1.007 can be excluded as it would 
imply T dust =13K. This identifies the line as CO(3-2) at 2=2.010. 

SPT0457-49: There is no line detected in the spectrum. The dust temperature would be « 22 K if the source is in 
the 2=1.74-2.00 redshift desert. As with SPT0128-51, a higher redshift with weak lines cannot be ruled out. This 
source would lie at 2=3.3 were its Tdust the same as the median of the unambiguous sample. 

SPT0459-58: A single CO line is detected at 98.40 GHz. If the line is identified as CO(4-3) at 2=3.685, the 
Ci( 3 Pi-> 3 P ) transition is in the band as well at 105.12 GHz . In this case the C i( 3 Pi ^ 3 P )/CO(4-3) flux density ratio 
limit is <0.15 (3er), comparable to the limit we observe for SPT0345-47. Therefore CO(4-3) cannot be excluded but 
would require an unusually low (but not unprecedented) C i/CO line ratio. CO(2-l) at 2=1.343 can be excluded based 
on the dust temperature (Td U st=14K). CO(3-2) at 2=2.514 implies Td U st=22K. The most plausible identification is 
CO(5-4) at 2=4.856 with T du st=41K. 

SPT0512-59: A single CO line is detected at 106.94GHz. CO(4-3) and CO(5-4) can be excluded as C i( 3 Pi -*. 3 P ) 
should have been detected given the bright CO line. CO(2-l) at 2=1.156 is unlikely as it implies Td us t=20K, but 
cannot be ruled out. Our preferred identification is CO(3-2) at z = 2.234 with Td U st=33K. 

SPT0550-53: A single bright CO line is identified at 111.67GHz. CO(2-l) at 2=1.064 is excluded (T du st=14K); 
for CO(5-4) at 2=4.160 CO(4-3) should have been detected. CO(3-2) at 2=2.096 and CO(4-3) at 2=3.128 are both 
plausible identifications with Td U st=22 and 31 K, respectively. 

SUPPLEMENTARY FAR-INFRARED PHOTOMETRY 

In this appendix, we show the supplementary FIR through mm photometric measurements used to determine dust 
temperatures, assign probabilistic redshift estimates to the sources with single-line detections and to show the repre- 
sentativeness of the dust colors of this subsample for the larger sample of 1.4 mm selected SPT sources meeting the 
same selection criteria. 

We used the LABOCA instruments at APEX to obtain 870 imaging. The observations took place during ESO 
and MPIfR o bserving time bet ween 2010 September and 2012 May. The observing strategy and data processing are 
described in (|Greve et al.ll2012l ). 

Herschel-SPIRE maps at 250, 350, and 500 /im were observed as part of program OT2_jvicira_5. The SPIRE data 
consists of a triple repetition map, with coverage complete to a radius of 5 arcmin from the nominal SPT position. 
The maps were produced via the standard reduction pipeline HIPE v9.0, the SPIRE Photometer Interactive Analysis 
(SPIA) package vl.7, and the calibration product v8.1. Photometry was extracted by fitting a gaussian profile to the 
SPIRE counterpart of the SPT detection and the noise was estimated by taking the RMS in the central 5 arcmin of 
the map which is then added in quadrature tothe absolute calibration uncertainty. 

For SED fits, we have added in quadrature an absolute calibration uncertainty of 10% for SPIRE, 15% for LABOCA, 
10% for SPT, and 10% for ALMA. 
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TABLE 4 
Far-Infrared and mm Photometry 
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Note. — Fluxes are given in units of mjy and include absolute calibration uncertainties. 2 mm & 1.4 mm fluxes have been deboosted. 
All other flux densities are photometric measurements at the ALMA position of the 1.4mm source. We note that source blending is 
typically not a problem for the photometry as the strong galaxy-galaxy lensing implies that the FIR light is dominated by a single lensed 
background object . Contamination by the lensing foreground galaxy can be ruled out by our ALMA high angular resolution 870/^m imaging 
HVieira et al.ll2013T ). 




Fig. 16. — Left: 870/^m LABOCA flux density as a function of 1.4mm SPT flux density for the 26 sources discussed in this paper (red) 
compared to the full sample of SPT sources which have been observed with LABOCA and Herschel-SPTRE (black). Right: Same as to the 
left but with 350/^m Herschel-SPTRE flux density as a function of 870/xm LABOCA flux density. 



